Unraveling the acoustic electron-phonon interaction in graphene 
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Using a first-principles approach we calculate the acoustic electron-phonon couplings in graphene 
for the transverse (TA) and longitudinal (LA) acoustic phonons. Analytic forms of the coupling 
matrix elements valid in the long-wavelength limit are found to give an almost quantitative descrip- 
tion of the first-principles based matrix elements even at shorter wavelengths. Using the analytic 
forms of the coupling matrix elements, we study the acoustic phonon-limited carrier mobility for 
temperatures — 200 K and high carrier densities of 10 12 — 10 13 cm -2 . We find that the intrin- 
sic effective acoustic deformation potential of graphene is H e g = 6.8 eV and that the temperature 
dependence of the mobility p ~ T~ a increases beyond an a — 4 dependence even in the absence 
of screening when the full coupling matrix elements are considered. The large disagreement be- 
tween our calculated deformation potential and those extracted from experimental measurements 
(18 — 29 eV) indicates that additional or modified acoustic phonon-scattering mechanisms are at 
play in experimental situations. 

PACS numbers: 81.05.Hd, 72.10.-d, 72.20.-i, 72.80.Jc 



I. INTRODUCTION 

Since the experimental realization of graphene,— its 
electronic properties and their understanding have been 
studied extensively both experimentally and theoreti- 
cally. 2- — While the intrinsic carrier mobility of graphene 
is predicted to be exceptionally high, the experimen- 
tal reality in substrate supported graphene involving 
charged impurities, electron-hole puddles, surface-optical 
phonons of the substrate, and disorder typically results 
in strongly reduced mobilities compared to the expected 
intrinsic value££ Together with scattering on acoustic 
phonons which manifests itself in a linear temperature 
dependence of the mobility at higher temperatures, these 
extrinsic scattering mechanisms typically dominate the 
mobility. The linear temperature dependence character- 
istic for acoustic phonon scattering has so far been ob- 
served in both supported^ - — and suspended 11 graphene 
samples. 

With the recent improvements in sample fabrication, 
the relative role of acoustic phonon scattering must 
be expected to become increasingly important in fu- 
ture devices. For example, samples with the commonly 
used SiC>2 substrate replaced by hexagonal boron ni- 
tride (h-BN) which has a lattice constant very close to 
that of graphene and an almost atomically flat surface 
with strongly reduced disorder ; 12 : 13 have shown highly 
improved transport characteristics with mobilities ap- 
proaching that of suspended graphene4 i 14 : 15 Further- 
more, the high energy of the surface-optical phonons 
of h-BN results in a significant reduction of surface- 
optical phonon scattering^ - — that for commonly used 
gate oxides starts to dominate the mobility around T ~ 
150 - 200 K££ 

When the mobility is dominated by acoustic phonon 
scattering, two transport regimes separated by the Bloch- 
Griineisen (BG) temperature Tbg = ZhkpCph/k^, can 



be identified^ Here, kp is the Fermi wave vector, 
c p h the sound velocity and fee the Boltzmann constant 
(Tbg ~ 57 K ^/n for the LA phonon with n measured 
in units of 10 12 cm -2 ). Since the BG temperature cor- 
responds to the acoustic phonon energy hujq — hcp^q 
for full backscattering at the Fermi level, short wave- 
length acoustic phonons are frozen out at temperatures 
T < Tbg restricting scattering processes to small scat- 
tering angles. The restricted phase space available for 
phonon scattering at T < Tbg results in a transition 
from the linear p ~ T behavior of the resistivity in the 
high-temperature regime to a stronger p ~ T a tempera- 
ture dependence in the BG regime where a = 4 (a = 6) 
in the absence (presence) of screening by the carriers 
themselves i 20 i 21 The BG behavior in the temperature de- 
pendence of the mobility (resistivity) has recently been 
observed experimentally X 

Existing theoretica l 20 ' 21 and experimental 7 - - —^ stud- 
ies of acoustic phonon-limited transport in graphene of- 
ten parametrize the interaction with acoustic phonons 
in terms of a coupling to a single effective acoustic 
phonon. The associated deformation potential coupling 
constant extracted from the experimentally measured 
temperature dependence of the resistivity range from 
- 18-29 eV&t&±L On the other hand, theoretical stud- 
ies of the acoustic electron-phonon coupling yield much 
lower values on the order of 3 — 4.5 eVi 22 ' 23 At the same 
time, different forms of the coupling matrix element are 
used in theoretical studie s 20 ' 22 making a direct compar- 
ison of the different values of the deformation potential 
difficult. 

Even though the effect of acoustic phonon scattering on 
the carrier mobility in graphene has been studied widely 
in the literature ) 20 ' 23- — a complete study considering the 
full details of the coupling matrix element is still lack- 
ing. The purpose of the present study is to provide a 
detailed analysis of the acoustic electron-phonon inter- 
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action in graphene and to establish the intrinsic value 
of the effective deformation potential. We shall focus 
on supported graphene where the flexural phonons are 
quenched and hence include only the transverse (TA) 
and longitudinal (LA) acoustic phonons. We use a first- 
principles method to calculate the electron-phonon inter- 
action^^ supported by the group-theoretical consider- 
ations of Refl30l. We then study the intrinsic phonon- 
limited mobility using a Boltzmann equation approach 
in the temperature regime — 200 K and for high carrier 



densities 



10 1 



10 cm where screening by the 



carriers suppresses other scattering mechanisms. Using 
the full coupling matrix elements, we find that a temper- 
ature dependence with a > 4 occurs even in the absence 
of screening. 



II. THEORY 

In the following, the carriers in graphene are described 
by massless Dirac fermions with linear dispersion e k = 
fwpk where vf ~ 1-0 x 10 6 m/s is the Fermi velocity. 
Within the Boltzmann equation approach j 6 ' 25 i 31 the mo- 
bility in graphene in the presence of (quasi) elastic scat- 
tering mechanisms is given by 
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where a xx is the conductivity, n is the two-dimensional 
carrier density and the density-of-states averaged relax- 
ation time (in units of time per energy) is defined by 



Tk- 



(2) 



Here, p(eu) — {g s g v /2-Kh 2 )e-^/ vp is the density of states 
of the graphene layer and g s — 2 and g v = 2 are the spin 
and valley degeneracies, respectively. At low tempera- 
tures and high carrier densities where ej? » /cbT this 
yields fx xx « ev 2 F T kF /e F . 

In the case of acoustic phonon scattering which can be 
treated as a quasielastic process, the relaxation time for 
each of the acoustic phonons (TA and LA) is given by 2 - 
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where #k,k' is the scattering angle and / k = /(£k) the 
Fermi function. The transition matrix element is given 
by 
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Here, gfc is the electron-phonon coupling, A denotes 
the acoustic phonon branches, and ftw q A the acoustic 



phonon energy. In the following, these quantities are cal- 
culated from first-principles. The phonons are assumed 
to be in equilibrium and populated according to the Bose- 
Einstcin distribution function Nq\ — N(hujq\). As scat- 
tering on both the TA and LA phonon is considered here, 
the total relaxation time for the K, If'- valleys is given 
by the sum of the individual phonon contributions as 
t -1 = J2 X rC" . As we show in the following section, the 
matrix elements of the electron-phonon coupling differ in 
the K and K' valleys (see e.g. Fig. [TJ. The Boltzmann 
equation must therefore be solved explicitly in both val- 
leys. In the absence of intervalley scattering which cou- 
ples the distribution functions in the two valleys, this 
can be done by considering the two valleys separately. 
In this case, the relaxation time entering the expression 
for the mobility in Eq. (flj becomes the valley-averaged 
relaxation time r = (r K + t k ,)/2, where t k/k , is the to- 
tal relaxation time in the individual valleys. Screening 
of the electron-phonon interaction by the carriers them- 
selve o 32 ' 33 has been considered elsewhere 2 ^ and will here 
be neglected. 

While analytic considerations have been given in 
RefsHalll we will in the present work resort to a nu- 
merical evaluatio n 34 ' 35 of the relaxation time in Eqs. Q 
and ([J|. This allows us to study the acoustic phonon- 
limited mobility in graphene with the full coupling ma- 
trix elements which have a more complex angular depen- 
dence than most often assumed (see e.g. Ref.20). The 
numerical approach also allows for a unified treatment of 
the high-temperature (T > Tbg) and Bloch-Gruneisen 
(T < Tbg) regimes. We note, however, that in the low- 
temperature regime where hujq ~ k^T, it is crucial that 
the phonon energy is retained in the Fermi function in 
Eq. ([3]) . In the high-temperature regime this requirement 
can be relaxed and the relaxation time can be put on a 
simple analytic form. 



III. INTERACTION WITH ACOUSTIC 
PHONONS 

In the following we use a first-principles DFT approach 
based on a fully microscopic description of both the elec- 
tronic Bloch states and the phonons to calculate the 
acoustic electron-phonon couplings in graphene i 28 ' 29 i 36 
The calculated values for the sound velocities c\ for the 
TA and LA phonons are reported in Table U together 
with other parameters used in this work. Due to their 
high phonon energies (> 100 meV), acoustic intervalley 
and optical phonons do not play a role in the considered 
temperature range and can therefore be neglected. 

The interaction between charge carriers and the acous- 
tic phonons in graphene can be written in the general 
form 



5kc 
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kq y 2Apujq X kq ' 
where A is the area of the graphene layer, p is the mass 
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density and Mj^ q = (k + q|(5VqA|k) is the coupling matrix 
element for scattering between the two Bloch states k and 
k + q due to a phonon with wave vector q, branch index 
A and frequency ujq\ = c\q where c\ is the sound velocity 
for the acoustic branches. The coupling is mediated by 
the change <5VqA in the microscopic crystal potential due 
to a unit displacement of the atoms along the mass-scaled 
normal mode vector e q >. Due to the full microscopic 
treatment of both electrons and phonons, Umklapp pro- 
cesses involving reciprocal lattice vectors are included in 
the coupling matrix element M£_. 

Figure Q] shows the absolute value of the coupling 
matrix elements M£ to the TA and LA phonons in 
the K, K' valleys as a function of the two-dimensional 
phonon wave vector q. The matrix elements in the two 
valleys are related through time-reversal symmetry as 
^kq A = M~k X ~q- This implies that carriers in the two 
valleys traveling in the same direction experience differ- 
ent electron-phonon couplings. It should be emphasized 
that the calculated matrix elements include the full sym- 
metry of both the electronic Bloch states and the phonon 
modes as given by the complete microscopic description. 

In order to emphasize the effect of the chirality of the 
carriers in graphene, the initial carrier state k is located 
on the right side of the Dirac cones 300 meV above the 
K, X'-points as indicated by the sketch in the top of the 
figure. As is evident from the figure, both the TA and 
LA phonons couple to the carriers with similar coupling 
strengths. While backscattering is suppressed for the LA 
mode, the situation is reversed for the TA mode where 
forward scattering is suppressed. In addition to suppres- 
sion of forward and backscattering, other directions with 
complete suppression of scattering also appears. This is 
a consequence of the inclusion of the symmetry of both 
phonons and electronic states. 



Parameter 


Symbol 


Value 


Lattice constant 


a 


2.46 A (LDA) 


Ion mass density 


P 


7.6 x 1(T 8 g/cm 2 


Fermi velocity 


Vf 


1.0 x 10 6 m/s 


Transverse sound velocity 


cta 


14.1 x 10 3 m/s 


Longitudinal sound velocity 


cla 


21.2 x 10 3 m/s 


Electron-phonon couplings 






Transverse 


Pta 


2.8 eV 


Longitudinal 


Qla 


2.8 eV 


Longitudinal 


/3la 


2.5 eV 


Effective coupling parameters 






Sound velocity 


Ceff 


20.0 x 10 3 m/s 


Deformation potential 




6.8 eV 



TABLE I: Material parameters for graphene used in the 
present work. The phonon related parameters have been ob- 
tained from first-principles as described in the text. The cal- 
culated sound velocities are in excellent agreement with the 
values reported in Refl27l. 
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FIG. 1: (Color online) Electron-phonon couplings to the 
acoustic TA and LA phonons in the K, if'-valleys (left and 
right columns, respectively) of graphene. The contour plots 
show the absolute value of the coupling matrix elements j M^ q | 
for a carrier energy of 300 meV as a function of the two- 
dimensional phonon wave vector q. The (white) circles corre- 
spond to k + q vectors lying on the constant energy surfaces 
of the Dirac cones given by the energy £k of the initial state 
in k as sketched in the top row. 



In the following, the first-principles coupling matrix el- 
ements are analyzed using the group-theoretical anal ysis 
of the electron-phonon interaction presented in Ref|30j. 
In the long- wavelength limit, the TA and LA phonons 
are strictly transverse and longitudinal, respectively^ 
and the electron-phonon interaction has a simple ana- 
lytic representation in the two-dimensional pseudospin 
basis^ Using the results of Refl30l. the coupling matrix 
elements can expressed in terms of the angles 6*k, # q and 
$k+q of the involved wave vectors. Including contribu- 
tions of order O(q), we find that the coupling matrix 
elements in the long-wavelength limit takes the following 
form in the if-valley, 



IM, 



TA I 
kq I 



qfi x 



sin 26>, 



(6) 



and 



M k L q A | = q x 



a cos 



'k+q 



+/3 cos 29, 



?k+q 



(7) 



for the TA and LA phonons, respectively. For the 
LA phonon, the first and second terms originate from 
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the deformation potential and the gauge field cou- 
pling mechanisms, respectively. The TA phonon cou- 
ples only through the latter . 3 °] 38 i 39 Since the interaction 
is Coulombic in nature, the overall couplings given in 
Eqs. © and ([7]) are here referred to as deformation po- 
tential couplings. 

With the coupling parameters listed in Table HI we 
find that the analytic expressions for the electron-phonon 
interaction in Eqs. © and to a high degree repro- 
duce the first-principles matrix elements for electron en- 
ergies up to ~ 750 meV. As the analytic coupling matrix 
elements are based on the phonon modes in the long- 
wavelength limit, the agreement is slightly worsened at 
shorter wavelengths where the mode vectors deviate from 
the long-wavelength modest This is most pronounced 
for the TA phonon. As the LA phonon retains its long- 
wavelength character far out in the Brillouin zone, the 
agreement between the coupling matrix elements here re- 
mains quantitative even at shorter wavelengths. In the 
BG regime where short wavelength phonons are frozen 
out, we note that it is the long- wavelength limit of the 
coupling matrix elements that governs the scattering of 
carriers. 

Often scattering on acoustic phonons is described by 
coupling to a single effective phonon mode with a cou- 
pling matrix element given b y 20 ' 25 

A<=~ cff(Z cos(^f±^ (8) 

where S e fj is the effective deformation potential and 
the angular part corresponds to the bare spinor overlap 
(Xk+q|Xk) of the electronic wave function. In contrast 
to the more complex angular dependence of the coupling 
matrix element predicted by the full microscopic treat- 
ment presented here, the angular dependence of the ef- 
fective coupling matrix element above suppresses only 
backscattering. In the high-temperature regime where 
equipartitioning of the acoustic phonons 7V q ~ k^T /Hlj^ 
applies, the relaxation time and the resistivity take the 
following simples forms 2 ^ 

1 1 . p = i3ML (9) 

r k h 3 Apv 2 F c 2 k ' P Ae 2 hpv 2 F c 2 ' { 1 

where the factor of 4 in the denominators stems from the 
chiral nature of the carriers through the assumed form 
of the coupling matrix element in Eq. ((HJ. These ex- 
pressions are used almost exclusively to extract the value 
of the effective acoustic deformation potential in experi- 
mental situations £i2rjL±l. 



IV. RESULTS 

In the following, we study the intrinsic acoustic 
phonon-limited mobility of graphene in both the BG and 
linear resistivity regime using the full coupling matrix 
elements as given by Eqs. (J5J) and ©. This allows us 




FIG. 2: (Color online) Inverse relaxation time (valley- 
averaged) for acoustic phonon scattering on the TA and LA 
phonon in the BG regime at n = 10 12 cm -2 (Tbg ~ 57 K). 
The full lines show the results obtained with the full matrix 
elements given in Eqs. ([6]) and Q and coupling constants ex- 
tracted from ab-initio calculations. The dashed lines show the 
result obtained with scattering on a single effective acoustic 
phonon with the coupling matrix element given by Eq. |(8j, 
a deformation potential of H e ff = 6.8 eV and sound velocity 
c = 20 x 10 3 m/s. 



to establish the value of the intrinsic effective acoustic 
deformation potential in graphene. 

In Fig. [2] we show the inverse of the valley-averaged 
relaxation time as a function of energy for different tem- 
peratures and a carrier density of n = 10 12 cm~ 2 corre- 
sponding to Ef ~ 117 meV and Tbg ~ 57 K for the LA 
phonon. Above the BG temperature, the inverse relax- 
ation time has the linear energy dependence of Eq. (|9]) 
and a slope proportional to the temperature. As the tem- 
perature is decreased below Tbg, the freezing out of short 
wavelength phonons and the sharpening of the Fermi sur- 
face result in limited phase space for phonon scattering 
and an increased life time of the carriers at the Fermi en- 
ergy. In the expression for the relaxation time in Eq. ([3]) 
this effect is accounted for by the Fermi and Bose dis- 
tribution functions. The limited phase space available 
for phonon scattering, manifests itself in the character- 
istic dip at the Fermi energy that evolves in the inverse 
relaxation time with decreasing temperature. 20 For all 
temperatures, the linear energy dependence of the high- 
temperature result in Eq. © is recovered in the e — > 
limit. 

By inspecting the individual contributions, we find 
that the inverse relaxation time to a large extent is dom- 
inated by the contribution from the TA phonon both in 
the high-temperature and the BG regime. The domi- 
nation of the TA phonon can be attributed to a num- 
ber of factors. From Eq. ([9]) it follows directly that the 
lower sound velocity of the TA phonon leads to higher 
scattering rate. Also related to the sound velocity is 
the lower BG temperature of TA phonon which allows 
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FIG. 3: (Color online) Mobility vs temperature for carrier 
densities 10 12 - 10 13 cm 2 (upper plot). The lower plot shows 
the temperature dependence of the exponent a in the temper- 
ature dependence fi ~ T~ a of the mobility. 



for full backscattering below the BG temperature of the 
LA phonon. Secondly, the coupling matrix element for 
the TA phonon allows for backscattering which is sup- 
pressed for the LA phonon. In the BG regime, the dom- 
ination of the TA phonon stems from the suppression of 
the coupling matrix element for the LA phonon in the 
long- wavelength limit (see Fig. [TJ . The observed domi- 
nance of the TA phonon is in contrast to the often used 
assumption that only the LA phonon couples to charge 
carriers in graphene^ 

In order to determine the intrinsic value of the effec- 
tive deformation potential in graphene, we also calcu- 
late the relaxation time using coupling matrix element 
in Eq. ([5]) . The dashed lines in Fig. [5] show the inverse 
relaxation time calculated with an effective deformation 
potential and sound velocity of 6.8 eV and 20.0 x 10 3 m/s, 
respectively. It is seen to reproduce the relaxation time 
based on the full matrix elements very well for the energy 
range shown. While the extracted value for the acoustic 
deformation potential is much smaller than experimen- 
tal values^^i it is in better agreement with recently 
reported ab-initio results yielding 4.5 eV^ 3 - 

Figure [3] summarizes the calculated acoustic phonon- 
limited mobility as a function of temperature for car- 
rier densities 10 12 — 10 13 cm -2 . The mobility calculated 
with the above-mentioned effective coupling parameters 



(dashed lines) reproduces the full calculation extremely 
well. For all carrier densities the mobility shows a transi- 
tion from the linear /i ~ T -1 high-temperature behavior 
to a more pronounced \x ~ T~ a temperature dependence 
with a > 1 in the BG regime. At temperatures T > Tbg, 
the decrease in the mobility with increasing carrier den- 
sity stems from the linear energy dependence of the den- 
sity of states which provides more phase space for phonon 
scattering for higher Fermi levels. At low temperatures 
T < Tbg where scattering on the full Fermi surface is 
frozen out, the local value of the coupling matrix element 
becomes important and the mobilities for the different 
carrier densities approach a common value. 

The right plot shows the temperature dependence of 
the exponent a for the same set of carrier densities which 
correspond to BG temperatures Tbg ~ 57 — 180 K for 
the LA phonon. From this plot it is more clear that 
the departure away from the linear temperature depen- 
dence happens for T ~ Tbg- For T < Tbg, the exponent 
increases monotonically. Surprisingly, the exponents ob- 
tained from the mobility calculated with the full coupling 
matrix elements do not saturate at a = 4 as predicted by 
the effective coupling matrix element (dashed lines) i 20 i 21 
Thus, even in the absence of screening, the mobility of 
graphene should take on a temperature dependence with 
a > 4 at sufficiently low temperatures. With carrier 
screening taken into account, this behavior is reinforced 21 

The purely intrinsic mobilities calculated here are sig- 
nificantly higher than previously reported theoretical val- 
ues.— This is reflected directly in the extracted defor- 
mation potential parameter which is considerably lower 
than commonly used values. For a carrier density of 
n = 10 12 cm~ 2 a room-temperature mobility in excess of 
10 5 cm 2 V -1 s _1 is here predicted. The associated scat- 
tering rate which is given by the relaxation time in Eq. Q 
in the high-temperature regime, is from Fig. [2] estimated 
to be r" 1 ~ 10 11 s _1 at the Fermi level corresponding to 
a mean-free path of A = vf /t ~ 1000 nm. Such an ex- 
tremely large mean-free path may open the opportunity 
to study coherent transport in relatively large graphene 
structures. We note that the quasiparticle scattering rate 
observable in photoemission spectroscopy (ARPES) de- 
velops a similar dip at the Fermi level for T < Tbg re- 
sulting in long-lived quasiparticles. 



V. CONCLUSIONS AND DISCUSSIONS 

In the present study the acoustic electron-phonon in- 
teraction in graphene has been analyzed in detail. The 
exact analytic forms of the coupling matrix elements 
in long-wavelength limit were found match the calcu- 
lated first-principles matrix elements almost quantita- 
tively even at shorter wavelengths. As previously pre- 
dicted,—^ the calculated mobilities showed a transition 
from a /j ~ T _1 to a /i ~ y— a temperature dependence 
with a > 1 below the BG temperature. However, con- 
trary to earlier studies we found that the full coupling 
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matrix elements cause the temperature dependence of 
the mobility to increase beyond a — 4 which is otherwise 
only observed when screening is included.— 

By fitting the results with an effective acoustic phonon 
the acoustic deformation potential is found to be 5 g = 
6.8 eV. Since this is much lower than the experimen- 
tally determined values 18 — 29 eV, our results suggest 
that the acoustic phonon-limited transport in substrate- 
supported graphene is at present not fully understood. 
Possible explanations for the large experimental deforma- 
tion potentials could be (i) substrate-induced modifica- 
tions of the band structur o 40 i 41 that modifies the chirality 
(and thereby the angular dependence of the coupling ma- 
trix element) of the electronic states and/or the Fermi 
velocity, and (ii) the existence of additional acoustic 
phonons not considered in the present work. This could 
be either the intrinsic flexural phonon which may be mod- 
ified when graphene is lying on a substrate or surface- 
acoustic phonons on the substrate which have previously 



been studied in 2DEGs42 The relatively large variations 
in the experimental deformation potentials also indicate 
that the deformation potential depends on experimental 
factors such as e.g. the substrate. 
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